Efficient Preparation of Quantum States With Exponential Precision 



Peter JaksclE 
Department of Physics and Measurement Technology, 
Linkóping Universïty, S-58183 Linkòping, Sweden. 
(Dated: February 1, 2008) 

It has been shown that, starting from the state |0), in the general case, an arbitrary quantum 
state |V>) cannot be prepared with exponential precision in polynomial time. However, we show 
that for the important special case when \ip) represents discrete vàlues of some real, continuous 
function ip{x), efficient preparation is possible by applying the eigenvalue estimation algorithm to a 
Hamiltonian which has ip(x) as an eigenstate. We construct the required Hamiltonian explicitly and 
present an iterative algorithm for removing unwanted superpositions from the output state in order 
to reach \tjj) within exponential accuracy. The method works under very general conditions and 
can be used to provide the quantum simulation algorithm with very accurate and general starting 
^S" states. 
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The first step of a quantum algorithm typically consists of resetting all registers to |0) followed by the preparation 
of an initial state \ip). In the general case, it has been shown impossible to prepare an arbitrary state with exponential 
precision in polynomial time 0. For special cases, however, tiris is still possible, e.g., when represents sampled 
oints of some continuous function ij)(x), an d certain integrals over tp(x) can be evaluated with exponential precision 

A useful application of such methods is to prepare initial states for quantum simulation Q. 
In this paper we show that quantum states that represent sampled, continuous functions can, under very general 
OO i conditions, be prepared efficiently. We also provide a method for doing this, using the quantum eigenvalue estimation 
algorithm by Abrams and Lloyd and the initial state preparation algorithm by Jaksch and Papageorgiou 0. 

Suppose that \if>) represents discrete vàlues of some real function i^{x). The intuitive idea is that if is an 
eigenvector of the (discretized) Schròdinger equation then it can be obtained by running the eigenvalue estimation 
l/"") ■ algorithm with an initial approximation generated by the technique in . The one-dimensional continuous Schròdinger 
equation for a particle with mass m, in a potential V(x), is defined as 
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4>{x)=e<f>{x), (1) 



where e is the energy. For convenience, we will choose units such that h 2 /(2m) = 1. By defining the potential as 

V[X) ~iP(x) dx 2 ' () 



ip(x) will apparently become an eigenvector with eigenvalue 0. Thus, if is measured from the eigenvalue estimation 
algorithm we know that the remaining (unmeasured) qubits will be in the state \ip). More precisely, assume that the 
continuous Hamiltonian (the expression in brackets in QJ) is discretized into a 2" x 2™ matrix H, and that U = e lHt 
can be implemented exactly on a quantum computer. An arbitrary eigenvector of U satisfics 

e^%) = e"'*^,) = e**^,). (3) 

Hence, A;i = 2-Ktpi + n27r, with n S Z and ipi € [0, 1). We want to prepare the specific state IV'o) for which Ao = 
(and ipo = 0). 

The discrete Hamiltonian H can be written H = T + V, where T is the second derivative operator, and V is 
the potential. The norm of the former is ||T|| = 2 2ïl 5j. We make the assumption that ||V|| = 0(2 2ïl ), so that 
< ||T|| + ||VIJ< 2 2n+p for some integer p, independent of n. For simplicity, we also assume that all elements 
of V are positive [8|, making H positively semidefinite, and that min^g A/ = 0(1), independent of n. By choosing 
t = l/2 2n+p all eigenvalues of U will apparently fall in the first quadrant, i.e., ifi < 1/4 for all l. 
Now, suppose that the eigenvalue estimation algorithm is started in the state 

2 n -l 

|0>|vI/ >= 53*10)1^), (4) 
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where the first register consists of only a single qubit. At the final step of the algorithm, before measurement, the 
state will be 

2"-l 

£ dí(flfo,0)|0>+ fl ( W ,l)|l>)M>. (5) 
Here, the function g is defined as (see 4J) 

{ sin(7i·(2y,-j))e' ri(y '- j72) 9 / ■ 
2sm(7r( V! -i/2)) ' Z W ^ J (g) 
1, 2<Pl=3- 

We remark that for the case j — 0, (fi ^ 0, g reduces to g((fii,j) = cosÍTrip^e^^ 1 . The probability of measuring the 
first qubit in the state |0), ensuring that \ïpo) is still in the superposition, is 

2™-l 

Pi=J2 \di\ 2 \g(w,o)\ 2 . (7) 

1=0 

Assuming that is measured, the second register will collapse to 

2"-l 

l*i) = 7p= E dig(<pt,0M)- (8) 
From induction it follows that repeated use of this procedure leads to 

p k = = \ £ K| 2 | 5 (^,o)r (9) 

= 7 =L= E *( 5 (^,o)) fc i^). (io) 

V i ' ' ' fe ;=0 

The probability of measuring every time (probability of success) is therefore 

2"-l 

Pf-P k =J2 \di\ 2 \g(w,0)\ 2 > Mo| 2 | 3 ((A),0)| 2 = Mo| 2 . (11) 
;=o 

The last equality follows since (^o = 0. With the method for initial state preparation in Q we can make |c£o| 2 > 1/2. 
From the expression for l^k) it is clear that only states with \g(<pi, 0)| close to 1 (ipt < 1/y/k) will survive as k grows. 
If we choose, e.g., k = 2n all states with ipi > 1/4 will have their amplitudes decreased by a factor more than 1/2™. 
By choosing a new value of t, t' = 2t we can reduce a second set of amplitudes exponentially. After 2n+p such steps 
we are left with 

P í 2n+P) = p( o) p( o) L·m p (2 K+p) El^| 2 |g^^)··-g(2 2 " + ^,0)| 2fc (12) 

I*i 2n+P) > = r= = = ; ; = E 0) • - ■ flf(2 2 "+P W , 0))*|^> (13) 

^ 1 (»)...p(»)...p 1 (^)...pf«+ri S 

2"-l 

P(success) = £ |di| 2 b(w,0)···5(2 2B +P W ,0)| 2 *> Mo| 2 |5(^o, 0) ■ • • 5 (2 2 "+^ , 0)| 2fc = |d | 2 . (14) 

1=0 

The condition that min/^o M — 0(1) ensures that at least one of the terms g(2^ipi,0) will satisfy gfòipijO) < 1/A, 
for some A > 1, independent of n and k. Hence, if k is polynomial in n the only term in the superposition that is not 
exponentially reduced is |^>o)- The algorithm can be implemented as a sequence of (2n + p)k circuits of the type in 
Fig. UI where = exp(iií2 s í). 
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FIG. 1: Schematic description of one step of the quantum state preparation algorithm. Here, Hd means the Hadamard gate; 
= exp(iH2 a t). When q reaches some value k it is reset to and s is increased with 1. 

The requirement that U can be implemented exactly can be weakened to requiring only exponential precision. In 
Fig. [l]the state before measurement is 



i|0)(J + ^W)|*W) + i|l)(/-l/W)|*W). 



The probability of measuring is Pq = \\I + C/^'|| 2 /4. For this outcome the second register will collapse to 



J=(J + l/W)|*W). 
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Suppose now that Ü^ 8 ' — + A can be implemented with exponential precision, i.e., ||A|| = 0(l/2 n ). For this 
case, we derive a lower bound on the probability of success. First, we make the following observation: 



p 9 (s) - v p [ q s) 



|(j + ^W)|*W)||-||(j + trW)|*W) 
< 1|| (Í/ W _ frW)|*W)|| = i||A|*( s ))|| < M. 

From JTJI it follows directly that |P ç (s) - P ç (;5) | < ||A||. Using this result we get 

p(o) . ,.p fc (2n+ P ) > (p fl» _ || A ||). . . (p^+p) _ ||A||) = PÍ 0) ■ ..p[ 2n+p \í - |0(||A||)|) & ( 2 "+^, 

which is bounded by a constant for exponentially small ||A|| [jj. Using H17|l and the definition of Pq^ we may also 
derive an upper bound on the difference between l'ï'ç+i) and l^^i) 
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1 "A|l< JèL = 0(||A||). 
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It can be shown ||| that the total error of the algorithm is at most the sum of the errors of the individual blocks. 
Hence, with the previous choice of k, 



I A tot || < (2n + p)feO(||A||) = 0(n 2 /2 n ) 



(20) 



We now consider the implementation of U. Suppose that we use an m:th order splitting formula (see 0) to 
approximate 
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where Wj are weight factors and Gj is either T or V. Defining the product on the right hand side LA 8 ' it follows from 
the definition of the matrix exponential that A is the Taylor expansion of U^ s ' — beyond m:th order. Provided 
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and e lTt can be 



that V(x) can be efficiently computed to exponential precision on a classical computer, both 
efnciently calculated on a quantum computer, as described in 0. 

The norm of A is of order 0((2 s t) m (\\T\\ + || V~||)™ 1 ). For some J3* > 1 we impose the (not very restrictive) condition 
\\T\\* + \\V\\* < -B*A„, where the norm is defined on a subspace spanned by a number of basis vectors \ipi) of H, 
containing not only the ground state, and A* is the maximum eigenvalue of these vectors ^(j- Defining B — sup-B*, 
where the supremum is taken over all subspaces, we choose t as t = l/(B2 n+p ). For k = — \n/ log 2 (7r/(4S))] all states 
with (pi > 1/(45) will have their amplitudes decreased by a factor more than 1/2™. 

We are now in position to modify the previous error analysis. If 2 s tB\* < 1/C, for some C > 1 independent of m, 
the error in f2*ï"Jl will be exponentially small, of order <3(1/C m ), on the subspace where A* is defined. 
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J2 d l (g( ( p l ,0)---g(2 s i p h 0)) k \^ 



2 s íBA,>l 



(22) 



For the eigenvectors in the sum in Q22)). if 2 s íA/ > 2ir, then at least one of the arguments 2^(pi will be O(l), and, hence, 
the amplitude of this state will be exponentially reduced. If 2 s tXi < 2tt, on the other hand, then 2 s ipi > l/(2wB). We 
make the assumption that B is bounded from above by a polynomial in n. Hence, 
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where D is greater than 1 independent of k. The fact that ||A|| = O(l) follows since ||C/ ( ' S ' ) || = ||£^ s - ) | = 1- 
choosing k as above, and m polynomial in n the right hand side of (|23|) can be made exponentially small in n. 



(23) 
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the same way it is possible to prové that ) — l 1 !'^)!! will also be exponentially small in n, and, thus, the same 



holds for the total error. In practice, one would have to guess the value of B. However, as long as the guessed value 
is larger than the true B it does not have to be very aceurate. 

For the sake of clarity, we have not yet discussed the issue of boundary conditions, but rather assumed implicitly 
that V'(O) = ^(1) = 0. This will be the case for most problems of practical interest and is also a requirement for the 
method in Q to work. Here, however, we will briefly outline a scheme for generalizing our method to functions with 
non-zero boundary conditions. One way of doing this is to extend the function to be prepared by smoothly attaching 
exponentially decaying tails on both sides. Practically, this can be achieved by adding, say, two extra qubits to the 
register holding the function to be prepared. As an example, states starting with qubits 00 will represent the left hand 
side tail; states starting with qubits 01 will represent the original function -0; and states starting with qubits 10 and 
11 will represent the right hand side tail. At the end of the preparation we will be left with the extended function. 
To get rid of the tails we can now measure the two qubits that we added to the register. With high probability this 
measurement will yield 01, implying that the extended state has collapsed to the original function ip. 

I am grateful to Anargyros Papageorgiou and Jan-Àke Larsson for very fruitful discussions and helpful comments. 



* Electronic address: petja@ifm.liu.se 
[1] E. Knill (1995), quant-ph/9508006. 

[2] C. Zalka, in Proc. Royal Soc. London (1998), vol. A 454, pp. 313-322, quant-ph/9603026. 

[3] D. S. Abrams and S. Lloyd, Phys. Rev. Lett. 83, 5162 (1999). 

[4] P. Jaksch and A. Papageorgiou, Phys. Rev. Lett. 91, 257902 (2003). 

[5] J. W. Demmel, Applied numerical linear àlgebra (Society for Industrial and Applied Mathematics, Philadelphia, 1997). 
[6] E. Bernstein and U. Vazirani, SIAM J. Còmput. 26, 1411 (1997). 
[7] A. J. Dragt, Phys. Rev. Lett. 75, 1946 (1995). 

[8] It is straightforward to relax this condition to allow for negative vàlues of V. 

[9] In fact, this is true even for polynomially small ||A|| 
[10] The inequality ||T|| + \\V\\ < \\H\\ +2||V[| is trivial. Suppose that ip is some eigenvector of H with eigenvalue A. For the case 
V > X, ip decays exponentially as e x ~ v . Hence, if A is the maximum eigenvalue on some subspace spanned by eigenvectors 
of H, the norm ||V|| will normally be of order A. The right hand side of the inequality will then also be of order A on this 
subspace. 



